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Abstract 


In  recent  years  the  lift-and-project  approach  has  been  used  successfully  within  a  branch- 
and-cut  framework  to  solve  large,  difficult  pure  and  mixed  0-1  programs  that  have  resisted 
solution  efforts  by  pure  branch  and  bound  codes.  The  approach  uses  a  linear  description  in 
a  higher  dimensional  space  of  the  convex  hull  of  the  di.sjunctive  set  created  by  imposing  oik* 
or  several  0-1  conditions.  By  solving  a  linear  program  derived  from  this  higher  dimensional 
representation  -  the  cut  generating  linear  program  (CGLP)  -  the  standard  lift-and-project 
procedure  obtains  a  deepest  cut  in  a  well  defined  sense.  We  propose  a  modification  of 
CGLP  that  allows  us  to  generate  not  just  one  deepest  cut.  but  a  class  of  cuts  with  desirable 
properties,  each  at  the  cost  of  one  extra  pivot  in  the  optimal  tableau  of  the  modified  CGLP. 


This  paper  investigates  a  certain  enhancement  of  the  lift-and-project  approach.  The 
underlying  research  has  its  immediate  background  and  starting  point  in  the  work  reported 
in  [5]  and  [6],  and  is  part  of  a  larger  project  pursued  jointly  with  Sebastian  Ceria  and  Gerard 
Cornuejols.  It  has  some  connections  with  the  work  of  Lovasz  and  Schrijver  [11]  on  matrix 
cones,  the  work  of  Sherali  and  Adams  [12],  but  primarily  with  the  work  on  disjunctive 
programming  [1,  2,  3]  done  in  the  70’s.  An  early  version  of  the  material  discussed  here  was 
circulated  under  [4].  For  some  related  work  on  lift  and  project  see  [7,  9,  10,  13]. 

The  lift-and-project  approach  to  mixed  0-1  programming  relies  primarily  on  the  following 
two  ideas  (results),  both  of  which  have  their  roots  in  the  work  on  disjunctive  programming, 
i.e.  on  optimization  over  unions  of  polyhedra,  done  in  the  seventies  ([!];  see  also  [3]  for  1 
below,  and  [2]  for  2  below). 


1.  There  is  a  compact  representation  of  the  convex  hull  of  a  union  of  polyhedra.  Namely, 
given  polyhedra  Pi  :=  {.t  G  iR"  :  A^x  >  6’}  ^  0,  i  G  Q,  the  closed  convex  hull  of 
Ui^qPi  is  the  set  of  those  x  G  J?”  for  which  there  exists  vectors  (y*,  j/q)  ^  i  E  Q 

such  that 


^  -i^Q)  =  0 

Ay-b%  >0  ieQ 

Po  >  0 

Eivo-^^Q)  =  1- 


The  number  of  variables  and  constraints  in  this  representation  is  linear  in  the  number 
[(^1  of  polyhedra  in  the  union. 


2.  The  closed  convex  hull  of  the  feasible  0-1  points  of  a  mixed  0-1  program 


min{c.T  :  Ax  >  b,  x  >  0,  x  .  e{0,i},  j  eNiCN} 


(1) 
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can  be  generated  by  imposing  the  0-1  conditions  successively.  Namely,  let  S\. . . .  ,Sp 

be  an  arbitrary  partition  of  A^.  and  for  ?  =  1 . p.  let  Q,  index  the  collection  of  0-1 

vectors  with  components  in  S,  (i.e..  |Q,|  =  2''“’’^).  Further,  let 

A'o  :=  {.r  G  1?"  :  A.r  >  b.  .r  >  0.  .7^  <1.  j  G  A’,} 
and  define  recursively 

Ki  conv{Ki--[  n  {.r  G  JR'"  ;  .r_,  G  {0. 1}.  j  G  5,}, 

7  =  1, . . .  ,p,  where  for  any  set  11’,  co7?r(ir)  denotes  the  clo.sed  convex  hull  of  IF.  Then 

Kp  =  covv{Ko  n  {.T  G  i??”  :  Tj  G  {0. 1},  j  G  A'l}). 

Conceptually,  the  lift-and-project  approach  is  supposed  to  work  as  follows.  Let  Ax  >  b 
denote  the  system  in  J?”  defining  A'o- 

0.1.  Solve  the  LP  over  A'o 

0.2.  Choose  5i  C  Ni  and  lift  the  problem  into  replacing  Ax  >  b  with  a  linear 

system  in  x  and  defining  the  convex  hull  of  A'o  H  {.r  :  Xj  G  {0. 1}.  j  G  Si}. 

0.3.  Project  the  system  onto  iF?"  (the  .r-space).  Let  o^.t  >  /?^,  h  G  T\.  be  the  .set  of 
projected  inequalities  that  together  with  Ax  >  b  define  A'l. 

Iterate  the  procedure.  At  step  t.  do 

t.l.  Solve  the  LP  over  A'f_i. 

t.2.  Choo.se  St  C  A'l  \u'l}S,  and  lift  the  problem  into  IP’^AQMn+^)  ^y  replacing  the  con¬ 
straints  of  Kt-\  with  the  linear  system  in  x  and  {y’.t/o},€Q,  defining  the  clo.sed  convex 
hull  of  A^,_i  n  {.T  :  Xj  G  {0. 1),  j  G  S,}. 
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t.3.  Project  the  system  onto  iR"  (the  .r-space).  Let  q^.t  >  h  G  Tf,  be  the  set  of 
projected  inequalities  that  together  with  those  of  Kt-i  define  Kf. 

After  p  iterations,  this  procedure  should  in  principle  yield  the  optimum  over  the  closed 
convex  hull  of  Kp,  i.e.  the  solution  to  the  mixed  integer  program. 

In  practice,  such  an  approach  is  not  workable,  since  at  every  iteration  the  number  of 
projected  inequalities  a^x  >  0^,  h  €  Tt,  is  exponential  in  n.  Therefore  the  standard  lift-and- 
project  procedure  of  [5],  [6]  generates  just  one  member  of  the  family  indexed  by  T/,  namely 
the  inequality  that  provides  a  “deepest  cut”  in  that  it  cuts  off  the  optimal  LP  solution  x 
over  Kt-i  by  a  maximum  amount. 

This  procedure,  suitably  embedded  into  a  branch-and-cut  framework  where  cuts  are 
generated  at  various  nodes  of  the  search  tree  and  made  globally  valid  through  a  cut-lifting 
step,  has  been  implemented  and  extensively  tested;  primarily  in  [6],  but  also  in  [7,  8,  14] 
among  others.  In  [6],  the  branch-and-cut  code  MIPO  using  lift-and-project  cuts  was  tested 
on  a  battery  of  29,  mostly  difl&cult,  mixed  0-1  programs  obtained  from  MIPLIB  and  from  the 
literature,  and  shown  to  be  able  to  solve  all  29  problems  in  computing  times  that  compare 
favorably  with  OSL,  CPLEX  and  MINTO.  In  [7]  the  same  approach  with  some  added  features 
was  tested  on  maximum  clique  problems  from  the  Second  DIMACS  Challenge.  In  [14]. 
S.  Thienel  compares  the  performance  of  his  ABACUS  branch-and-cut  code  in  two  different 
modes  of  operation,  one  using  lift-and-project  cuts  and  the  other  using  Gomory  cuts,  with 
the  outcome  that  the  version  using  lift-and-project  cuts  is  considerably  faster  on  all  hard 
problems,  where  hard  means  reqiiiring  10  or  more  minutes.  Finally,  in  [8]  the  authors  report 
that  in  the  framework  of  a  parallel  branch  and  cut  code  which  uses  a  variety  of  cuts,  adding 
lift-and-project  cuts  for  the  harder  problems  helped  solve  them  faster. 

As  mentioned  earlier,  the  standard  lift-and-project  approach  generates  one  deepest  cut 
whenever  it  solves  a  cut-generating  linear  program.  The  objective  of  the  research  discussed 
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here  is  to  find  a  middle  ground  between  the  two  extremes  of  generating  all  the  inequalities 
Q*.T  >  h  G  Tt.  and  generating  just  one  deepest  cut.  This  middle  ground  is  defined  as 
generating  all  those  inequalities  of  the  family  indexed  by  T,  that  are  tight  for  the  optimal 
solution  X  of  Kf. 

Consider  the  mixed  0-1  program  (1)  and  its  linear  programming  relaxation 

min{c.r  :  Ax  >  b.  x,  <1,  j  E  Ni,  x,  >0.  j  E  N} 

r  -  '  ‘  (LP) 

=  min{c.r  :  Ax  >  b}. 

Let  S  C  Ni  and  define 

Ps  :=  conv  {.r  E  IR"  :  Ax  >  b.  Xj  E  {0. 1}  for  j  E  S}. 

Then  there  exist  inequalities  a^x  >  ,  h  E  T,  such  that 

Ps  :=  {.r  €  1?"  ;  a^r  >  h  E  T}, 
but  these  inequalities  are  exponentially  many  and  not  readily  available. 

Problem.  Given  x  such  that 


cx  =  min  {c.7-  :  x  E  Ps}. 

find  the  inequalities  q^x  >  0^.  h  eT  CT  that  are  tight  for  x. 

Application.  If  S  is  ‘small",  x  can  be  found,  for  instance,  by  branch  and  bound.  Solving 
the  above  Problem  would  then  enable  us  to  replace  (LP)  by  min  {cx  :  Ax  >  b,  a^x  > 
0^ ,  h  E  T'}.  whose  solution  x  has  Xj  E  {0,1},  j  E  S.  and  repeat  the  procedure  for  some 
S"  C  Ni\  S.  In  other  words,  solving  the  above  Problem  is  a  way  of  ‘  consolidating"  the 
results  of  a  partial  branch  and  bound  run  into  a  tighter  LP  relaxation. 

Discussion.  Let  Q  index  the  collection  of  all  0-1  vectors  with  components  in  5,  i.e.  |Q|  = 
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Consider  the  extended  formulation  of  Ps  in  n  +  |(5|(^  +  1  )-dimensional  space: 

Ps  =  {.T  G  iR”  :  X  -  =  0 

ie<5 

Ay*  -  'b*yo  >  0,  ieQ 
E?/o  =  1  }• 

ieQ 

Here  6*  :=  where  6^  is  defined  as  follows.  If  the  original  right  hand  side  is  written  as 
b  :=  ^^*0)5  where  6^  €  {0,  —  with  the  first  n  components  equal  to  0  and  the  remaining 

p  components  equal  to  —1,  then  6'  corresponding  to  the  t-th  0-1  vector  indexed  by  Q  is 
obtained  from  6^  by  changing  the  —1  in  5*^  to  0  in  6^  for  all  those  variables  that  have  to 
be  forced  to  0,  and  changing  the  0  in  6°  to  1  in  6*  for  all  those  variables  to  be  forced  to  1. 
Notice  that  the  constraints  Pq  >  0,  i  ^  Q,  have  been  omitted,  as  they  are  implied  by  the 
other  inequalities. 

To  express  P5  in  .r-space,  consider  the  projection  cone 


f 

a  —  u'A 

=  0 

W  :=  < 

{a,  {u*}iGQ,  P) 

u^b’  — 

>  0 

p  =  0  i  e  Q  > 

The  ineqiialities  q^.t  >  defining  Ps  in  .r-space  correspond  to  the  (a,/3)-components 
of  the  extreme  rays  of  W.  The  task  is  to  identify  those  extreme  rays  of  W  that  give  rise  to 
inequalities  a^x  >  (3^  satisfied  at  equality  by  x. 

Solution.  Define  the  linear  program 

max  f3  —  ax 

—  a  -f  u'A  =  0 

P  -  =  0  ieQ  (LDl(.r)) 

<  k 
»eQ 

tt’  >  0,  ieQ 

where  k  >  0  and  e  =  (1, . . . ,  1). 


5 


Theorem  1  The.  inequalities  o^.r  >  0^  of  P5  satisfied  at  equality  by  r  are  precisely  those 
corresponding  to  the  [a,  0)- components  of  optimal  solutions  to  (LDl{r))  such  that  o  ^  0. 

Proof.  The  feasible  set  of  (LDl(.f))  is  the  cone  W  truncated  by  the  inequality  <  k. 

a  normalization  device.  Basic  solutions  to  (LDl(.f))  .such  that  o  0  are  in  1-1  correspon¬ 
dence  with  extreme  rays  of  IT.  Consider  now  the  linear  program  dual  to  (LDl(.f)). 

min  IryH 

-'Ey’  = 

ieQ 

Ay'  -  Vyl  -f  ei/o  >0.  i  (LPl(x)) 

Evo  =  1 

»ec? 

yfl  >  0. 

(LPl(.r))  has  an  optimal  solution  of  the  form  y'‘  =  .t.  y'^^  =  1  for  the  particular  €  Q  that 

corresponds  to  .f  (i.e.  for  the  unique  such  that  D'’t  >  d'"),  and  {ff  .ylf)  =  (0.0)  for  all 
i  E  Q  \  yo  =  0.  Hence  (LPl(.f))  has  an  optimal  solution  with  value  0.  Further,  notice 
that  can  be  pivoted  into  the  basis  with  value  0.  i.e.  without  changing  the  solution.  It 
follows  that  (LDl(.f))  has  an  optimal  solution  of  value  0  for  with  =  k.  i.e.  (LDl(.r)) 

has  a  nonzero  optimal  solution  with  value  0.  The  (o.  /?)-component  of  such  a  solution  defines 
an  inequality  satisfied  at  equality  by  .r  if  and  only  if  o  ^  O.D 

Thus  one  way  of  solving  our  problem  is  to  generate  nonzero  basic  optimal  solutions  to 
(LDl(.r))  satisfying  =  k.  Although  (LDl(.f))  is  large,  it  has  a  strong  structure  that 

can  be  exploited. 

Working  in  a  subspace.  We  should  mention  at  this  point  that,  as  in  the  case  of  the  lift- 
and-project  procedure  de.scribed  in  [5].  (LDl(.r))  can  be  re.stricted  to  the  relevant  subspace. 
and  the  inequalities  obtained  in  this  fashion  can  be  lifted  into  the  full  space  by  the  same 
technique  as  in  [5].  For  our  purposes  F,  the  index  set  of  the  relevant  subspace,  will  be  defined 
slightly  differently,  so  as  to  always  contain  S.  although  the  components  of  .r  indexed  by  S 
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are  not  fractional.  Also,  we  will  assume  w.l.o.g.  that  for  j  €.  S.  Xj  =  1  implies  that  xj  is 
basic;  and  for  j  ^  S,  if  Xj  is  nonbasic  then  xj  =  0  (i.e.  all  nonbasic  Xj  at  their  upper  bound, 
if  any,  have  been  complemented).  We  define 

F  :=S[J{j  eN\S  :  xj  >  0}, 

and  we  will  denote  by  xp  the  subvector  of  x  with  components  indexed  by  F.  The  matrix 
Af  is  obtained  from  the  matrix  A  by  deleting  the  columns  indexed  hy  N  \  F  and  the  rows 
corresponding  to  the  inequalities  Xj  >  0  and  —Xj  >  —1  for  j  E  N\F.  Further,  the  vectors  bp 
and  bp,  i  €  Q,  are  obtained  from  b  and  6’,  i  £  Q,  by  removing  the  components  corresponding 
to  the  inequalities  Xj  >  0  and  -Xj  >  -1  for  j  E  N  \  F.  The  problem  (LDl(S))  can  then  be 


replaced  by 

max  /?  — 

axp 

(X  +  Ajr  = 

0 

P 

-  u'b'jp  = 

yu‘e  < 

i^Q 

u’  >0,  i  E  Q 

0  i  E  Q 

k 

(LDI(xf)) 

where  a  and  i  E  Q,  have  been  redimensioned  according  to  Ap  and  are  to  be  read  as 
shorthand  for  ap  and  Up,  i  E  Q. 


Elimination  of  (a, /?).  Now  let  €  Q  be  the  index  associated  with  x,  i.e.  the  (unique)  i 
such  that  Ay'  >  b'  is  satisfied  by  x.  Using  a  =  u'’Ap  and  /3  =  ^’*6^  eliminate  a  and  /?, 

and  imposing  equality  in  the  last  inequality,  (LDI(.tf))  can  be  restated  as 
min  u'*  {ApXp  —  b'p ) 

u'' Ap  —  u' Ap  =  0 

-u'4’p  +  u%  =  0  (LDI(xf)) 

^u'e  =  k 
ieQ 

u'  >  Q,i  E  Q 

Proposition  2  A  feasible  solution  to  (LD1(.tf))  is  optimal  if  and  only  if 

{Apxp)j  >  {Fp)j  =»  uif  =  0. 
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Proof.  From  the  foregoing  discussion,  the  optimum  of  (LDl(.f/.-)j  has  value  O.D 


While  (LDl(.r;r))  has  fewer  variables  than  (LD1(.77  )).  this  latter  formulation  has  the 
advantage  of  having  the  values  of  o  and  d  readily  available  when  needed.  Further,  since 
q.I3  are  unconstrained  in  sign,  they  can  be  made  part  of  the  .starting  basis  when  .solving 
(LDl(.f/'))  and  kept  in  the  basis  till  the  end.  so  that  their  presence  does  not  affect  the 
number  of  pivots.  Therefore  we  will  prefer  to  work  with  (LD^.ff^)). 

Generating  inequalities.  It  is  ea.sy  to  see  that  optimal  solutions  to  (LDl(.r/. ))  are  highly 
dual  degenerate,  i.e.  there  are  typically  many  different  optimal  solutions.  In  fact,  from  the 
proof  of  the  above  Theorem,  it  follows  that  all  the  reduced  costs  corresponding  to  uj  for  all 
j  and  i  e  Q  \  {f.}  are  0,  since  these  reduced  costs  are  the  slack  variables  of  the  inequalities 
A’y'-b%>0oi  (LPl(.f)K). 

It  should  be  noted  at  this  point  that  while  the  complete  .set  of  inequalities  o'*.?-  > 
h  G  T' .  tight  for  .r.  obviously  cuts  off  the  optimal  solution  .f  to  (LP).  this  is  not  nece.ssarily 
the  case  for  each  individual  inequality  indexed  by  V.  In  order  to  generate  inequalities  that 
are  guaranteed  to  cut  off  .f.  one  may  add  to  the  constraints  of  (LDl(.T/r))  the  inequality 
0  —  axf  >  e  with  some  small  e  >  0.  One  advantage  of  doing  this  is  that  it  also  guarantees 
Q  0.  Thus  generating  the  inequalities  o.r  >  /?  of  P5  that  are  tight  for  .f  can  be  done  by 
the  following  procedure: 

1.  Solve  (LDl(.f|r)). 

2.  Find  alternate  optimal  solutions  by  pivoting  into  the  basis,  one  at  a  time,  nonbasic 
variables  with  zero  reduced  cost  for  which  the  minimum  ratio  (of  the  primal  simplex 
method)  is  nonzero. 

3.  For  each  new  vector  (q./5)  obtained  in  this  way.  use  the  ‘mininuim  required  angle” 
criterion  to  discard  ineqtialities  too  clo.se  to  some  earlier  ineqtiality. 
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4.  Lift  each  inequality  into  the  full  space. 


An  efficient  implementation  of  this  procedure  may  take  advantage  of  the  special  structure 
of  LDl(.Ti7).  Although  this  linear  program  is  relatively  large  ([(Jl  +  1  times  the  size  of  the 
original  LP  restricted  to  the  subspace  of  the  variables  indexed  by  F),  on  the  other  hand  is 
highly  structured.  Indeed,  the  coefficient  matrix  of  LD1(.tf)  has  the  general  form 

IQ  I  blocks 


0 

-I 

0  ••• 

0  ' 

0 

-I 

0 

••• 

0 

0 

-I 

0 

0  ••• 

Al 

1 

0 

0 

0 

1 

0 

0 

_52T 

0 

1 

0 

0 

0  ••• 

-IIQIT 

where  the  |Q|  blocks  in  the  upper  part  of  the  matrix  are  all  identical,  and  the  differences  are 
restricted  to  the  last  |Q|  rows  of  the  matrix,  containing  the  right  hand  side  of  each  term  in 
the  disjunction. 

Note  that  in  the  absence  of  any  special  structure  |Q|  =  i.e.  the  number  Q  of  terms  in 
the  disjunction  underlying  LDl(.rf')  is  exponential  in  the  number  l^l  of  0-1  variables  being 
arbitrated.  However,  the  presence  of  some  striicture  can  cut  down  this  number  drastically, 
in  some  cases  making  |Q|  linear,  rather  than  exponential,  in  |5|.  For  instance,  if  the  system 
At  >  b  contains  an  inequality  of  the  form  x{S)  <  1,  i.e.  a  generalized  upper  bounding 
constraint,  then  |Q|  =  |5|  -I-  1,  since  .ts,  the  subvector  of  x  with  components  indexed  by  S, 
is  either  one  of  the  |5|-dimensional  unit  vectors  or  the  1 51-dimensional  zero  vector. 
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Example  1. 


2.7-1  + 

•'?'2 

+ 

3.7-;, 

+ 

•7.1 

4.ri  + 

4.r2 

+ 

3.7;i 

— 

.7.1 

> 

3 

-.Ti  + 

2.7-2 

— 

2.7-;) 

+ 

3.7-4 

> 

1 

.Ti  - 

X2 

+ 

2.7-3 

+ 

to 

> 

1 

.T  >  0, 

(o  i).  i 

=  1, 

,2.3 

The  optimal  solution  to  the  LP  relaxation  (in  which  Tj  €  {0. 1}  is  replaced  by  0  <  .Tj  <  1, 
j  =  1,2,3)  is  X  =  (0.578.  0.270.  0,  0.346). 


We  choose  S  :=  {1,2}.  i.e.  impo.se  the  condition  Xj  €  {0.1}  for  j  =  1,2.  and  we  find 
the  optimal  solution  x  =  (0, 1,0, 1)  to  the  LP  relaxation  subject  to  this  condition.  Thus  we 
define 


F  :=  SU{jeN\S:  x^  >  0} 

=  (1,2.4), 


and  work  with  the  problem  (LDl(.rf )): 
min  (3  —  ax 


P 

P 

P 

P 


Q 

a 

a 

Q 


AptP 


+  AZtP 


+  ApiP 


+  Ap^P 


-  FpW 


-  bjTv^ 

-  bfiP 

-  bf'iP 

,iP.  iP,  u'*  >  0 


0 

0 

0 

0 

0 

0 

0 

0 
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Here 


Xp  =  (.r,..r2..r4)^  =  (0. 1. 1)^ 


\T 
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and 


4 

4 

-1 

3 

3 

3 

3 

-1 

2 

3 

1 

1 

1 

1 

1 

-1 

2 

1 

1 

1 

1 

Ap  = 

1 

1 

1 

6].= 

bl  = 

1 

1 

II 

-0 

1 

1 

-1 

-1 

-1 

-1 

-1 

-1 

where  the  blanks  represent  zeros. 

Solving  (LDl(.T/r))  yields  an  optimal  solution  whose  (a,/?)  components  are 

ai  =  0.559,  02  =  0.102,  04  =  1.119,  /3  =  1.221. 

The  lifting  coefficient  for  03  is  easily  seen  to  be  0,  and  dividing  through  with  f3  to  obtain  a 
right  hand  side  of  1  yields  the  cut 

0.458.T1  +  0.084.T2  +  O.9I6.T4  >  1. 


Performing  several  pivots  in  the  optimal  tableau  according  to  the  rules  discussed  above 


yields  the  additional  optimal  solutions 

£1  £2 

0.875  0.448 
0.654  0.218 
0.976  0.650 

In  each  case,  the  lifting  coefficient  for  Q3  is  0 
following  three  additional  cuts: 

O.623.T1  -|-  O.319.T2 
0.429.T1  +  O.I43.T2 
O.6OO.X1  +  0.400.r2 


Qu  0 
0.957  1.407 
1.307  1.525 
0.976  1.626 

and  dividing  through  with  /3  we  obtain  the 

+  0.68 1.T4  >  1 

+  0.857.r4  >  1 
+  O.6OO.T4  >  1 


Each  of  the  above  four  cuts  is  satisfied  at  equality  by  x.  Each  of  them  happens  to  cut 
off.f.D 
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Duplication  of  inequalities.  One  difficulty  with  the  procedure  outlined  above  is  that  it 
will  frequently  generate  multiples  of  the  same  inequality.  This  can  happen  in  two  different 
ways: 

1.  A  pivot  may  change  some  or  all  of  the  variables  t/’.  i  G  Q.  but  leave  the  components 
of  Q,/?  unchanged. 

2.  A  pivot  may  change  the  components  of  [a.  6)  but  change  them  proportionally,  i.e.  it 
may  leave  the  ratios  Oj//5.  j  G  N .  unchanged. 

Both  kinds  of  difficulties  can  be  handled  by  an  appropriate  check  of  the  pivot  column,  say  g. 
To  avoid  problem  1.  one  has  to  make  sure  that  g  has  a  nonzero  entry  in  at  least  one  of  the 
rows  associated  with  the  components  of  (o./?).  To  avoid  problem  2.  one  has  to  make  sure 

that  the  ratios  Oj/po, . fi/gn  are  not  all  equal.  Here  (gii)  denotes  the  entry  of 

g  in  the  row  containing  o,  (3). 

It  is  not  always  po.ssible  to  avoid  a  pivot  which  changes  the  n’  but  leaves  the  cut  o.r  >  /? 
es.sentially  unchanged.  If  necessary,  such  a  pivot  must  be  performed  in  order  to  obtain 
another  working  basis  that  yields  access  to  a  new  set  of  pivots. 

Generating  facets  of  Ps.  Under  what  circtimstances  do  the  inequalities  generated  by  the 
above  procedure  define  facets  of  P5? 

To  answer  this  question,  we  note  first  that  when  Ps  is  full  dimensional,  then  the  reverse 
polar  cone 

-  /5)  ^  :  ar  >  13  for  all  .r  G  Ps} 

of  Ps  is  pointed;  which  means  that  P|  is  the  conical  hull  of  its  extreme  rays.  Assume  for  the 
time  being  that  this  is  the  case.  Then  the  facets  of  P5  are  in  a  one  to  one  correspondence 

d 

with  the  extreme  rays  of  Pg. 
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On  the  other  hand,  from  the  representation  of  Ps  in  given  at  the  beginning 

il 

of  our  discussion,  it  follows  that  Pg  also  has  a  higher  dimensional  representation,  namely 


>  0, 

i  G  Q, 

such  that 

II 

(q,  /?)  € 

Q  — 

u'  A 

=  0 

13  - 

u'V 

=  0,  ieQ  ^ 

The  constraint  set  defining  this  higher  dimensional  representation  differs  from  that  of 
LD1(.t)  only  in  the  absence  of  the  normalization  inequality.  Adding  this  inequality  replaces 
the  cone  Pg  with  a  polyhedron 


Q  ” 

u'  A  =  0 

(a,;5)GiR”+' 

/5  - 

u‘0  =  0  i  E  Q 

E 

VI 

i£Q 

u'  >  0,  i  e  Q  J 

whose  extreme  points  (vertices)  are  in  one  to  one  correspondence  with  the  extreme  rays  of 
the  cone  Pg,  hence  with  the  facets  of  Ps- 

But  the  polyhedron  .P5  is  nothing  but  the  projection  onto  the  subspace  of  the  (a,/3) 
variables  of  the  feasible  set  of  (LDl(x)).  Thus  the  question  that  we  started  out  with  reduces 
to  the  following:  when  does  an  extreme  point  (basic  solution)  (q°, /3°,  i  E  Q)  of  the 

feasible  set  of  (LD1(.t))  project  into  an  extreme  point  (a°,/?°)  of  P5? 

One  answer  to  this  question  is  that  the  (a°,  /?°)-component  of  a  basic  feasible  solution  to 
(LD1(.t))  is  an  extreme  point  of  .P5  if  and  only  if  there  exists  a  linear  function  of  (q,  /?)  that 
attains  its  unique  maximum  at  (a,/?)  =  (a°, y0°).  A  basic  feasible  solution  to  (LDl(.f))  that 
satisfies  this  criterion  will  be  called  regular  (see  [1]  for  a  discussion).  To  check  the  regularity 
of  a  basic  feasible  solution  to  LDl(.f)),  one  may  want  to  carry  a  second  objective  function 
row.  in  which  the  vector  (1,  .f )  is  replaced  by  one,  if  it  exists,  that  makes  the  (q  ,  /?)  component 
of  an  optimal  solution  unique,  i.e.  produces  nonzero  reduced  costs  for  all  nonbasic  variables 
whose  pivoting  into  the  basis  would  affect  the  (a, /5) -component.  Methods  for  finding  such  a 
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vector  or  showing  that  none  exists  are  currently  under  investigation.  So  far  we  have  assumed 
that  Ps  is  full  dimensional.  When  this  is  not  the  case.  i.e.  dim  —  d  <  v .  then  is  not 
pointed,  hence  has  no  extreme  rays;  and  the  dimension  of  its  lineality  spaee  L  is  ??  —d.  In  this 
case  there  is  a  one  to  one  correspondence  between  facets  of  Ps  and  {n  —  d  +  l)-dimensional 
faces  of  Pg.  Since  the  latter  are  not  easy  to  handle,  it  is  preferable  in  this  case  to  work  with 
P5  n  Z/-*-,  where  is  the  orthogonal  complement  of  the  lineality  spaee  L  of  P|.  n  is 
a  pointed  cone  whose  extreme  rays  are  in  a  one  to  one  correspondence  with  the  {n  —  d  +  1)- 
dimensional  faces  of  P^,  hence  with  the  facets  of  P,.  The  analytical  expre.s.sion  for  this  cone 
is 


P|nL^  =  Pq,^)  G  E 


rj  +  l 


a.r  >  /?  for  all  .r  £  Pg 
Q7’  =  0  for  all  V  £  L 
where  “all  v  £  L"  can  be  replaced  by  “all  r  £  B{L)"  for  some  basis  B{L)  of  L.  Note  that 

dimZ/  =  n  —  d  and  thus  a  basis  B{L)  has  v  —  d  elements. 


When  it  comes  to  the  polyhedron  P|.  the  normalization  inequality  <  k  effec¬ 

tively  bounds  all  the  variables,  including  the  components  of  {a.0)\  and  thus  the  smallest 
dimemsional.  i.e.  {71  -  d  +  l)-dimensional.  faces  of  P|  corre.spond  to  extreme  points  of  P|. 
However,  this  correspondence  is  no  longer  one  to  one.  as  the  same  facet  of  Pg  can  now  be 
defined  by  different  inequalities  a.r  >  0,  corresponding  to  different  extreme  points  of  P^-. 
Nevertheless,  among  all  the  equivalent  inequalities  a.r  >  0  defining  the  same  facet  of  Pg. 
there  is  only  one  (modulo  a  multiplicative  factor)  whose  normal  lies  in  L-^ .  If  we  amend  the 

u 

constraint  set  of  P|  by  the  system  of  equations 

ar  =  0  for  all  r  £  B{L), 

i.e.  replace  P|  with  PIdL-^,  we  make  sure  that  the  correspondenee  between  extreme  points 
of  P5  n  L-*-  and  facets  of  Pg  are  one  to  one. 


We  illustrate  this  last  point  on  an  example. 


Example  2.  The  traveling  salesman  polytope  P  defined  on  the  complete  directed  graph  G 
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with  4  nodes  has  dimension  5  (=  number  of  variables  (12)  minus  rank  of  equality  system 
(2  X  4  -  1)).  The  6  (affinely  independent)  tours  of  G  are  (1,2, 3, 4),  (1,2,4, 3),  (1,3, 2, 4), 
(1, 3, 4, 2),  (1, 4, 2, 3),  (1, 4, 3, 2).  One  of  the  facets  of  P  is  defined  by  the  inequality 

Xi2  +  .Ti3  +  .^32  +  .T34  +  .T43  <  2, 

known  as  T2,  which  is  satisfied  at  equality  by  5  of  the  6  tours,  the  exception  being  (1, 4, 2, 3). 
However,  the  same  facet  is  also  defined  by  each  of  the  lifted  3-cycle  inequalities 

.Ti2  +  X24  -h  .'^41  +  2.T21  <  2, 

^13  +  •'^34  +  ^41  +  2.T43  <  2, 

^■24  +  ^43  +  ^32  +  2.T34  <  2, 

which  are  easily  seen  to  be  satisfied  at  equality  by  exactly  the  same  5  tours  as  the  T2 
inequality.  However,  none  of  the  above  four  inequalities  has  its  normal  vector  in  the  subspace 
L-*-  generated  by  the  tours.  Indeed,  every  vector  in  satisfies  all  the  eqiiations  of  the  form 

outdegree  of  z  =  outdegree  of  j 
indegree  of  z  =  indegree  of  j 
outdegree  of  z  =  indegree  of  j 
indegree  of  z  =  outdegree  of  j, 

just  as  every  tour  does:  the  normals  of  these  equations  are  all  in  L,  the  lineality  space  of  P^. 
However,  none  of  the  four  inequalities  defining  our  facet  satisfies  this  system  of  equations. 
The  unique  inequality  that  defines  the  same  facet  satisfies  the  system  is 

3.ri2  +  2xi3  -f  3.T21  +  2x24  +  2.T32  +  3x^4  -I-  2.T41  +  3.T43  <  8. 


Indeed,  it  is  easy  to  check  that  this  inequality  (i)  is  valid,  (ii)  is  satisfied  at  equality  by 
the  same  5  tours  as  the  above  four  inequalities,  and  (iii)  satisfies  the  system  of  equations 
defining 
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